Decoherence under many-body system-environment interactions: a stroboscopic 
representation based on a fictitiously homogenized interaction rate 
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An environment interacting with portions of a system leads to multiexponential interaction rates. 
Within the Keldysh formalism, we fictitiously homogenize the system-environment interaction yield- 
ing a uniform decay rate facilitating the evaluation of the propagators. Through an injection proce- 
dure we neutralize the fictitious interactions. This technique justifies a stroboscopic representation 
of the system-environment interaction which is useful for numerical implementation and converges 
to the natural continuous process. We apply this procedure to a fermionic two-level system and 
use the Jordan- Wigner transformation to solve a two-spin swapping gate in the presence of a spin 
environment. 
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I. INTRODUCTION 



The control of open quantum systems has a fundamen- 
tal relevance for fields ranging from quantum information 
processing (QIP) [J to nanotechnology 0,0, Q. Typi- 
cally, the system whose coherent dynamics one wants to 
manipulate, interacts with an environment that smoothly 
degrades its quantum dynamics. This process, called "de- 
coherence", can even be assisted by the own system's 
comp lexity [||. Since environment induced decoherence 
S 0, constitutes the main obstacle towards QIP, a 
precise understanding of its inner mechanisms [2, H, Gil 
is critical to develop strategies to control the quantum 
dynamics. 

The usual way to obtain the dynamics is to solve a 
generalized Liouville-von Neumann differential equation 
for the reduced density matrix. There the degrees of free- 
dom of the environment are traced out to yield a quantum 
master equation (QME) A less known alternative 

is provided by the Keldysh formalism [l2| in the inte- 
gral representation proposed by Danielewicz (l3j . On one 
side, it uses the well known perturbation to infinite order 
in selected terms provided by the Feynman diagrams. On 
the other, this integral representation has the advantage 
of being able to profit from a Wigner representation for 
the energy-time domain. This last representation is par- 
ticularly meaningful in the fermionic case since it allows 
ones to define energy states and their occupations simul- 
taneously with the physical time [l4j]- In that case, one 
can transform the Danielewicz equation into the general- 
ized Landauer-Buttiker equation (GLBE) [lj,[l5( to solve 
the quantum dynamics of the system. When the system- 
environment (SE) interaction is spatially homogeneous, 
i.e., it has an equal interaction with each component of 
the system, the dynamics becomes particularly simple 
because there is a uniform SE interaction rate. However, 
there are many situations where one should incorporate 
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multiple rates as different subsets of the system could 
suffer diverse interaction processes. While this might 
not possess a great challenge to the evaluation of steady 
state transport properties, in quantum dynamics, one is 
confronted with what appears to be a much more diffi- 
cult problem. Here, we present a procedure to convert 
a nonhomogeneous problem with multiple SE interaction 
rates, into one that has a common rate. Through a rein- 
jection procedure, we instantaneously neutralize the fic- 
titious decays restoring the populations and, eventually, 
the coherences. In order to illustrate the procedure, we 
apply this technique to a model that represents a single 
fcrmion that can jump between two states while an exter- 
nal fermionic reservoir is coupled to one of them. This 
provides decoherence due to a through space Coulomb 
interaction and can feed with an extra particle through 
tunneling processes. While the parameters and approxi- 
mations involved in this model are especially designed to 
be mapped to a problem of spin dynamics, it could also 
be adapted to represent a double quantum dot charge 
qubit [16j. In that case a double dot is operated in the 
gate voltage regime where there is a single electron which 
can jump between the two coupled states, where only 
one of these states is coupled to an electron reservoir. 
This inhomogeneous SE interaction yields a multiexpo- 
nential decay rate. We introduce fictitious interactions 
to obtain a common interaction rate which leads to a ho- 
mogeneous non-Hermitian effective Hamiltonian. In the 
specific model considered, we analyze how different SE in- 
teractions, e.g., tunneling to the leads and through space 
Coulomb interaction, modify the quantum evolution. A 
particular advantage of the fictitious symmetrization is 
that it leads naturally to a stroboscopic representation 
of the SE processes. This leads to a very efficient numer- 
ical algorithm where the quantum dynamics is obtained 
in a sequence of time steps. Finally, we resort to the 
Jordan- Wigner mapping between fermions and spins to 
apply the procedure to a spin system. This allows us to 
give a first-principle derivation of the self-energies used 
in the stroboscopic model introduced in Ref. [17| to ex- 
plain the puzzling experimental dynamics observed [l8j 
in a spin swapping gate [l9l ]. 
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II. SYSTEM 

Let us consider an electron in a two-state 
system asymmetrically coupled to an elec- 
tron reservoir, as shown in Fig. [Ha), with 
the total Hamiltonian TL = Tt s + H E + Hse- 
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FIG. 1: (Color online) (a) System-environment (SE) repre- 
sentation. Dashed circles and solid circles represent the sys- 
tem and the environment states respectively. Dashed lines 
are hopping interactions while wiggly lines are through-space 
Coulomb interactions, (b) Self-energy diagram summing up 
the different interactions with the environment in a local ba- 
sis. The lines with arrows are exact Green's functions in the 
absence of SE interactions. The double dashed line represents 
the effective SE interaction, (c) Retarded Green's function at 
site 1. The interaction with the environment is to infinite or- 
der in the self-energy given in (b). (d) Particle density func- 
tion at site 1. The double dashed lines represent the effective 
interactions local in time and space summed up to infinity 
order. 



The system Hamiltonian is 

Hs = E x cXc x + E 2 c\c 2 - V12 (c^c 2 + c^Ci) , (1) 

with c\ the standard fermionic creation (destruction) 
operators. The Ei are the energies of the i-th local state 



whose spin index is omitted. The hopping interaction 
V12 gives the natural frequency ujq = 2V\2/h of the tran- 
sition between the states 1 and 2. The environment has 
a similar Hamiltonian, 

00 00 
n E = ]T - Vij (c+S + cpi) , (2) 

i=3 M=3 

where the Vy determines the topology of the interac- 
tion network in the environment states. The system- 
environment interaction is described by 

+ ^23 ' ^2a^3a^ta^2a ~ ^23 iptcfisa + ^ta^2a) } ) (3) 

The first two terms on the rhs represent the Coulomb 
interaction of an electron in site 2 with an electron in 
site 3, the first site of the reservoir. U 2 ^ lr is the standard 

direct integral and U^^'^ is the small exchange integral 
which we include for completeness. The third term is the 
hopping interaction between sites 2 and 3. 

III. SYSTEM EVOLUTION 
A. Quantum dynamics in the Keldysh formalism 

We are interested in the study of the evolution of an 
initial local excitation in the system. Let us consider the 
initial excitation with a particle on site 2 and a hole in 
site 1 which is described by the non-equilibrium state, 

|* n . c .) = C+C, |tf eq .> , (4) 

where |\& C q.) is the thermodynamical many-body equilib- 
rium state at high temperatures which is the regime of 
NMR spin dynamics. In this condition I'J'oq.) is the mix- 
ture, with equal weight, of all the possible Slater deter- 
minants [20j . The evolution in a complete norm preserv- 
ing solution is described by the particle and hole density 
functions 

G< (*2,*l) = \ (*n.c.| ct (fi) C % (t 2 ) |* n . c .) (5) 

and 

Gfj (fc.tl) = ~ (*„.c.| C, (t 2 ) Cj (ti) |* n . e .) , (6) 

that describe spatial and temporal correlations. In these 
expressions, the creation and destruction operators are 
in the Heisenberg representation. Notice that in contrast 
with the equilibrium definitions of Gf^ (t2,ti), now they 
have an implicit dependence on the initial local excita- 
tion. The probability amplitude of finding a particle in 
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site i at time t 2 when it initially was in site j at time t% is 
described by the retarded Green's function of the whole 
system 

G R (ia.fi) =e{t 2 ,h) [G> (f 2) fi) - G< (t a ,t a )] 

= [G^tLfa)]*. (7) 

The reduced density function G < (f, f), where matrix in- 
dices are restricted to i,j <G {1,2}, is equivalent to the 
single particle 2x2 density matrix and G R (fa, f i) is an 
effective evolution operator [2l| ■ If the system is isolated, 
the Green's function in its energy representation is ob- 
tained by a Fourier transform (FT) with respect to the 
time interval 5f = fa — fx, 



G 0R (e, t) = J G 0R (f + §5f, t - \5t) exp[i£<5f//l]d5f, 
where f = ^ (fa + f i) . In a time independent system 



(8) 



G (e,f) = G 0K (e) = [eI-H s ]- 



(9) 



After including SE interactions, the Green's function 
defines the reduced effective Hamiltonian and the self- 
energy S R (e) [23, 



H 



eil' 



( e ) = sI-[G R (e)] _1 



He 



S R (e). (10) 



Here, the exact perturbed dynamics is contained in the 
nonlinear dependence of the self-energy S R on e. For 
infinite reservoirs the evolution with H g. is nonunitary, 
hence, the Green's function has poles at the "eigenener- 
gies", e„, that have imaginary components [23j], 

-2Im£ R (£„)/ft=l/T SE = 2r SE //i. (11) 

These account for the "decay rates" into collective SE 
eigenstates in agreement with a self-consistent Fermi 
golden rule (FGR) H]. Similarly, Re£ R (ej = Ree v - 
represent the "shifts" of the system's eigenenergies £°. 
The evolution of the density function for the reduced 
open system is described using the Keldysh formalism 
12 , [l3| . The density function in the Danielewicz form 
11 is 



G< (f 2 , fi) = fi. 2 G R (f 2 , 0) G< (0, 0) G A (0, h) 

dt k dUG R (f 2 , t k ) £< (f fc , *,) G A (t,, ti) . 

(12) 

The first term is the "coherent" evolution while the sec- 
ond term contains "incoherent reinjections" through the 
self-energy function S < . This compensates any leak from 
the coherent evolution reflected by the imaginary part of 
E R (see [0]])- We remark that this expression is valid for 
a noncorrelated initial state which is our case of interest. 



For a correlated state see Ref. [25[ . The key to solve Eq. 
(fT12j) is to build up an expression for the particle (hole) 
injection and retarded self-energies S < ' > J (fi,f2) and 

S R (f l7 < 2 ) = 0(f 1 ,f 2 )[S > (f 2 ,f 1 )-S < (f 2 ,f 1 )]. (13) 



For this purpose, we use a perturbative expansion on 
TfsE like that used in Ref. [26j for the Coulomb inter- 
action and in Ref. for the hopping interaction. The 
first order in the perturbation expansion is the standard 
Hartree-Fock energy correction which does not contribute 
to S < because it is real. We focus on the second-order 
terms, with Feynman diagrams sketched in Fig. [IJb). 
The injection self-energy is 



l^ 23 | 2 ^ 2 G| (t k , U) G< (t,, t fc ) Gf 2 (t k , ti) S i2 S 2j 

+ \V 23 \ 2 Gg(t k ,t l ) S i2 S 2j , (14) 



where U- 



23 



-2U. 



is the net Coulomb in- 



teraction between an electron in the system and one in 
the reservoir. The direct term contributes with a fermion 
loop and an extra spin summation which is represented 
in the —2 factor (l3| . The first term in Eq. (jT4j) corre- 
sponds to the direct and exchange self-energy diagrams 
shown in the last line of Fig. [ljb). The first two dia- 
grams schematize the creation of an electron-hole pair in 
the environment and its later destruction. The last term 
in Eq. (fT4"|) and the last diagram of the same figure is the 
hopping to site 3 which allows the electron to perform a 
full exploration inside the reservoir. To take into account 
the different time scales for the dynamics of excitations 
in the system and in the reservoir, we use the energy- 
time variables: the physical time t\ = \ (f& +fi)i an d 
the domain of quantum correlations 8t\ — t k — U. This 
last is related to an energy e through a FT Thus, 
in equilibrium, 



G 3 < 3 (£,f i ) = i27riV 3 (£) f 3 (Mi), 



C33 



(Mi 



-i27r7V 3 (£) [l-f 3 (M0] 



(15) 
(16) 



where N3 (e) is the local density of states (LDoS) at the 
surface of the reservoir. Assuming that the environment 
stays in the thermodynamical equilibrium and k^T is 
much higher than any energy scale in the bath (high tem- 
perature limit), the occupation factor is 

f 3 (Mi)=f 3 - (17) 
Fourier transforming on e one obtains 

G33 {ti + fA-f)= i2^ 33 (Sti) f 3 (18) 

and 

G&fa+f ,fi-f ) = -i27r a, (ft,) (l-f 3 ), (19) 



where 



de 



53(5*0= / JV 3 ( £ )e^-. 



(20) 
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Replacing in Eq. (M 



Sti leading to 



|[/ 23 | 2 ft 2 (2^) 2 [ ff3 (^)] 2 f 3 [l-f 3 ] 

X ^22 (*i + if 1 ! *i — <^2^2j 



±|V 23 | Z i27T 53 (5*0 ( iJfJ ( 21 ) 

where the (\i 3 f 3 ) associates f 3 with E< and (1 — f 3 ) with 
E>. 

In summary, we are left with the task to evaluate the 
time dependent self-energies and the integral in Eq. (fl"2"| . 
We will focus on the parametric regime corresponding to 
the experimental conditions of the spin swapping gate. 



B. Environment in the wide band or fast 
fluctuations regime 



As occurs with the generalized Landauer-Biittiker 
equations for linear transport, an essential ingredient is 
the possibility to assign a Markovian nature to the en- 
vironment. We are going to see that this appears nat- 
urally from the formalism when the dynamics of exci- 
tations within the environment is faster than the time 
scales relevant to the system. In order to separate the 
different physical time scales involved in the problem, we 
start changing to the energy-time variables in Eq. (fT2|) . 
Evaluating in t% = ti = t, the integrand becomes 



t r-t 

o J-t 

2 



x G« ( Mi + Sk.) £< + «i - f ) G A (t; - f ,t) 



(22) 



The environment unperturbed Green's function <? 3 (St{) 
decays within the time scale H/Vb where Vb is the charac- 
teristic interaction inside the reservoir. In the wide band 
regime (Vb S> V12) h/Vs becomes much shorter than 
the characteristic evolution time of G 22 (*i + U ~ 
given by fi/Vi 2 . Then, as the main contribution to the 
integral on Sti of Eq. (fT2|) is around the time scale H/Vb 
we can replace G^ 2 (ij + - by Gf 2 (ti,t{). Fol- 
lowing the same assumption we replace G R (t,ti + 
by G R (t, and G A (ti - 4fS t) by G A (t b t) . In this fast 
fluctuation regime, only Eg (tj + 4fs ti — ^) depends on 



(*i) 



E^ 



2 ' Ci 



2 , 



d5£i 



|C/ 23 | 2 n 2 (27T) 2 



[33 (<J*i)f dJti 



f 3 [1 - f 3 ] 



x G\\ (t h ti) 5 i2 5. 



2.; 



i27T 



ff3 (^i)dfti 



2j 
1-13 



(23) 



which is local in space and time. Here, because of the 
limit V12/VB — > 0, the correlation function of site 3 
becomes a representation of the Dirac delta function. 
Thus, any perturbation at site 3 is almost instanta- 
neously spread all over the environment (as compared 
with the time scale of the system dynamics) and hence 
the occupation at site 3 remains constant. This assump- 
tion for the time scales can be seen in Fig. [ljb) as a 
collapse of a pair of black dots, along a vertical line, into 
a single point. This justifies the expansion of Fig. [ljc) 
and the use of the ladder approximation containing only 
vertical interaction lines in Fig. |T](d) . 
A generalized decay rate is given by 

l/r SE (e, ti) = 2r SE (e, ti) /h = -2 Im E R (e, f;) /H (24) 

= i[E 2 A 2 ( e ,t i )-E« 2 (e,t i )] (25) 

= i[E> ( £ ,ti)-E< (eA)], (26) 

where 



Eg (MO 



E S (ti + f , ti - 5 -t) e iE5u/h dSti. (27) 



We start assuming Ei = for i = 1, ..,00, so the only 
relevant energy scale of the system is Vi 2 Vb- As men- 
tioned above, in the wide band limit the correlation func- 
tion gs (Sti) becomes a representation of the Dirac delta 
function. In this way, the term e leSt U h f the integrand 
of Eq. (|2"T)) is evaluated for St x — giving a value equal 
to 1. Thus, using Eq. ([23]) . we obtain for the decay rate 
which in the wide band limit is constant in time and in- 
dependent of energy, 



^WB l[ E 22(*i)- S 2 2 (*i)] 



(28) 



= \U. 



2:S 



I 1^1 



(Tt 



(27T) 2 

f r v ) 



[<? 3 (<5ti)] 2 d<5ti 



ff3 (Sti) dSti 



f 3 [1 - {3} 



(29) 
(30) 



where we have used t ^> H/Vb to equal E^(e,ti) = 



Eg (ti) and define 



h\u, 



2:; 



2vr 2 



[ 53 (Sti)} 2 dSti 



f 3 [1 - f 3 ] 



(31) 
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the Coulomb decay rate, and 



IV* 



93 (Sti) dSU 



(32) 



the hopping decay rate. If one assumes that the environ- 
ment ([2]) can be represented by a linear chain with near 
neighbor hoppings equal to Vb and Ei = 0, the LDoS is 
(seeRef. ^) 



N 3 (e) = 1/(ttV b ) 
Thus, the Green's function 
33 (SU 



e 
2Vr 



i himu 



2irV B 



St, 



(33) 



(34) 



is proportional to the first-order Bessel function and de- 
cays within a characteristic time H/Vb- Assuming that 
fs = 1/2 and the integration limits in the T's expressions 
are taken to infinity because t ~ h/V%2 ^> H/Vb (wide 
band approximation), one obtains 



frv-f P2 3 \ 2 ^r 



and 



|r v = ¥ IV 



2 1 

J;l ttVb' 



(35) 



(36) 



Since the interaction is local in time, the reduced density 
results as follows 



G< (i, t) = h 2 G R (i, 0) G< (0, 0) G A (0, t) 

+ f dt i G R (t,t i )-£ < (t i )G A (t i ,t) 
Jo 

which is complemented with 



£< ( ti ) 





2T u hG<(t u t i ) + 2T v h(j-i 3 ) 



(37) 



(38) 



Here, the propagators G R (t, 0) and G A (0, t) that enter 
in both terms are obtained from the effective Hamiltonian 
of the reduced system, 



H, 



-Via 
— V12 — irsE 



(39) 



where Tse is energy independent and the assumption 
Ei = for all i assures that the self-energies are purely 
imaginary. This effective Hamiltonian is obtained from 
Eq. (TIT))) and using Eq. (|13[) previously converted into 
the energy-time variables. This means that we first 
change the variables (iijia) to (ii,<5ii) and then based 
on a Fourier transform on Sti we obtain the energy-time 
representation which results in the independence on both 
e and U. 



The above procedure results in an equation of the form 
of the GLBE. However, the Hamiltonian is asymmetric 
in the SE interaction complicating the form of the as- 
sociated propagator. The apparent complexity to solve 
this equation contrasts with the homogeneous case where 
the evolution of the GLBE was obtained (l5( through a 
Laplace transform. Our strategy will be to induce such 
a form of symmetry. 



C. Fictitious homogeneous decay 

The main difficulty with Eq. (|37[) is that it involves 
multiple exponentials. In order to create propagators 
with an homogeneous decay, i.e., a single exponential 
factor, we introduce fictitious interactions Sg c with the 
environment. The symmetric Hamiltonian becomes 



H 



off. 



-Via 
-V12 — iLsE 

-i^rgE — V2 

— V12 — i^SE 



-Ur 



2-SE 





SE 



(40) 



Here includes the fictitious interactions which, in 

the present case, produce a leak of probability in site 1 at 
a rate Vs^/h while in site 2 inject probability at the same 
rate. Both states of H sym . interact with the environment 
independently with the same characteristic decay rate 
rgE/ft. Note that this rate is half the real value. The 
propagators of Eq. (fr2|) have now a simple dependence 
on t as 



G R (t, 0) = G 0R (t, 0) e 



where 



and 



G??(t, 0)=G° R M) = ~ h cos (^t) 



(41) 



(42) 



Gi R (t, 0) = G^it*) °)* = ^ sin ( y*) (43) 

are the isolated system propagators. The reduced density 
evolution is now, 

G< {t, t) = h 2 G 0R (t, 0) G< (0, 0) G 0A (0, t) e -*/(2r SE ) 

+ f'dtiG^^tOKym (ti)G OA (i i ,t) e -^ t 0/(2.s E ) j 

Jo 

(44) 

which is similar to the GLBE [lj, [l5| . It is easy to see 
that the introduction of negative (positive) imaginary 
parts in the diagonal energies of the effective Hamilto- 
nian produces decay (growth) rates of the elements of the 
density function which, being fictitious, must be compen- 
sated by a fictitious injection self-energy 

EfiU (*0 = ( s fic* + >"U) G< (U, U) . (45) 
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In our case, this results in an injection that includes the 
compensation effects for the symmetrized interaction 



£< ym , (tO = £< (ti) + S< . (tO 





o 2r^(if 3 ) + 2r £/ ^G< (t i5 to 
r SE ?iG< (t b tO o 







-r SE ?iG< (ti,to 



(46) 



Here, the second term is proportional to the system den- 
sity functions G^ (ti,t;) injecting and extracting density 
on sites 1 and 2, respectively, to restore the real occupa- 
tion. It is important to remark that the escape Ty given 
by V23 in Hamiltonian © or the process of current leads 
of Refs. [3, HH are only compensated at a constant rate 
by the reservoirs. In this case, the injection self-energy 
is proportional to the density function in the environ- 
ment. In contrast, for voltage probes, elect ron-phonon 
self-energies (as in Ref. [lj]) or our Coulombic Ty require 
an immediate charge compensation. The same is true for 
the fictitious processes and this is indeed the situation of 
Eq. (|46[) where the injection self-energy is proportional to 
the instantaneous system density function. Thus, the fic- 
titious injection self-energy compensates instantaneously 
the fictitious leak and injection of Eq. (|40|) . This is more 
easily seen once Eq. (|4"4")) is integrated into a Trotter-type 
form, as we will discuss in connection to Eq. (|76p . The 
symmetrization method is, in essence, a redistribution of 
terms in the evolution equation (|12|) that has a simpler 
resolution. 

We can rewrite the last expression to separate the pro- 
cesses that involve density relaxation (through injection 
and escape processes) and pure decoherence (through lo- 
cal energy fluctuations) as follows 



S< m , (ti) = £< (ti) + £< (ti) 



ir 



SE 



2py 











(f 3 -fG<(ti, tO) 

iG<(t u ti) 

fG<(ti,ti) 



Here 



and 



?G< (ti,*i) = ? / G< ( £ ,tO^=f 2 (tO 
1 1 / 2irn 



-G< (ti, to =f x (to 



(47) 



(48) 



(49) 



while, remembering that according to Eqs. (|24[) and ([50]) . 

Tse = IV + ry , we define 



Pv = Ty/rsi 



(50) 



as the weight of the tunneling rate relative to the total 
SE interaction rate. Since the initial state has the site 2 



occupied we have that 

jGfj (0,0) = S i2 S 2j . (51) 
1 J 

Introducing Eq. (|47|) into Eq. (|44|) and using 

— = f r SE , (52) 

TSE 

we get two coupled equations for GJ\ and G^ 2 as follows: 

T^ii (*>*) = |^G?^(t,0)| 2 e-* /(2TSE) + 
\hG™ (t, ti)| 2 e -(*-*0/(a^.) 2pv ^- [f 3 - f G< (UA)} 

+ J |fi G ?f (t,ti)| 2 e-^)/(^ E )^L[fG<(ti,ti)] 

+ htvG^itAfe-^^p- [fG< (MO] , 

(53) 



5G< (t,t) = |^(t,0)| 2 e-*/( 2 -) + 
^(t.tOl'e-C*"*')/^) 2p y [f 3 - fG< (ti,ti)] 

+ / |^f(t,ti)| 2 e-^)/(^)-^[fG<(ti,ti)] 

+ J |nG» (t,ti)| 2 e^)/^)^[fG< (ti,ti)]. 

(54) 

In each equation, the first term is the probability that a 
particle initially at site 2 be found in site 1 (or 2) at time 
t having survived the interactions with the environment 
with a probability e~ t ^ 2TSE K The second term describes 
the process of injection/escape of particles enabled by 
the hopping from/towards the reservoir, where the last 
of such processes occurred in the time range (ti, t; + dtj) 
with a probability 2pv ^^r- The injection (escape) at site 
2 fills (empties) the site to level it up to the occupation 
factor f 3 . The third and fourth terms take into account 
the last process of measurement at time t; due to the SE 
interaction with a probability 2^:- This confirms our 
interpretation that in Eq. (|47|) the dissipation processes 
are in (t) while (t) involves pure decoherence. It 
is clear that by iterating this formula, one gets a series 
in the form represented in Fig. [IJd) . 

In summary, Eqs. (|53|) and (T54|) are valid within the 
following assumptions: (a) The system is in the high tem- 
perature limit V±2, U23, V23, Vb «C k^T . (b) The environ- 
ment is assumed to have a very fast dynamic as com- 
pared with that of the system Vi2,£723,V23 <§; Vb (fast 
fluctuation regime). This is achieved through a specific 
model of the environment with a very wide band where 
this property shows up as a flat and broad local den- 
sity of state at the "surface" site 3. These two are the 
central assumptions. The consequences of them are the 
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following limits: (c) The self-energy is described by the 
self-consistent second-order term, (d) The environment 
remains in equilibrium, (e) The retarded and advanced 
Green's functions are calculated from a non-Hermitian ef- 
fective Hamiltonian which is independent of energy and 
constant in time. 

The previous conditions allow us to develop a different 
strategy for the solution of the spatially inhomogeneous 
evolution equation (GLBE): We fictitiously symmetrize 
the effective Hamiltonian to impose an homogeneous de- 
cay of the coherent dynamics. Consequently, we compen- 
sate the resulting artificial injections and/or leaks based 
on a fictitious part in the injection self-energy. 



D. Dynamics of a swapping gate 



is purely real or imaginary i.e. lot] = 0. Also, 



The solution of the coupled Eqs. ([55)1 and ([54")) involves 
a Laplace transform. We consider a parameter range 
compatible with the spin problem where f3 < 1 while 
we allow the tunneling relative weight py in the range 
[0, 1]. In a compact notation, the density function results 
as follows 



iG<(t,t) = l-a t 



- Rot - oi cos [(w + irf) t + tpo] e~ Rlt . 

(55) 

Here, the decay rates Rq, R%, and rj, and the oscillation 
frequency u> are real numbers associated with poles of 
the Laplace transform. The amplitude ao is also real 
while, when u> = 0, the amplitude a\ and the initial phase 
<£>o acquire an imaginary component that warrants a real 
density. These observables have expressions in terms of 
adimensional functions of the fundamental parameters in 
the model. Denoting 



x = w t se , 

and remembering that 

pv = rv/r SE , 

we define 

cj)(pv,x) = i (x 2 - p v - i (l-p v f j . 

and 



(56) 
(57) 

(58) 



X (pv,x) = {4(1- Pv ) (9x 2 -2(l-pv) 2 + 18 Pv ) 



+ 12 



(4x 6 -((l-pv) 2 + 12p v ) x 4 
+4p v (b{l-p v ? +3p 2 v) x 2 

1 s. 

-4p 2 v ((1 - p v f - P 2 V ) 



The observable "frequency" , 

fl , . 6(p v ,x)\ 



(59) 



(60) 



Rq 



. 4>(pv,x) 1 I \ 1 

>— 7 T - ~X(PV,X) +p V + o I 1 ~PV) , 

X(Pv,x) 6 3 J T SE 

(61) 

* = (62) 



and 



ao = 



1 2 (uj 2 - t? 2 ) + 2R 2 - ujj 
2^-^) + (R -R 1 ) 2 - 
1 



a 2 



«3 



2 (w + 17?) 

(2^^! - LP 2 ) (Rq - gl) + 2 (lj 2 - 7? 2 ) flp 

(c 2 -,y 2 ) + (i? - J R 1 ) 2 
1 c^g + 2i?§ - 
2( W 2 - 7? 2 ) + ( J R -i 1 > 1 ) 2 ' 



a i =a 2+ a 3' tan (0o ) 



02 
03 



(63) 

(64) 

(65) 
(66) 



The oscillation frequency u> in Eq. (|60|) has a critical 
point a; c at a finite value of x showing a quantum dy- 
namical phase transition for which u> and r\ in Eq. (|55[) 
exchange their roles as being zero and having a finite 
value, respectively. A full discussion of this issue for a 
spin system is presented in Ref. [17] . Here, the dynami- 
cal behavior changes from a swapping phase to an over- 
damped phase. This last regime can be associated with 
the quantum Zeno effect [28| where frequent projective 
measurements prevent the quantum evolution. Here, this 
is a dynamical effect [29T |30T | produced by interactions 
with the environment that freeze the system oscillation. 

Figure [5] shows typical curves of fG< (t,t) in the 
swapping phase. The different colors correspond to dif- 
ferent SE interactions, pv = 0, 0.5, and 1, which are 
Coulomb (IV = 0), isotropic (IV = IV) , and pure tun- 
neling (IV = 0) interactions rates. The hopping inter- 
action does not conserve the net energy in the system 
inducing a dissipation which is manifested through the 
nonconservation of the number of particles in the sys- 
tem. This is the case of pv 7^0 where the final state 
of the system has the occupation probability of the sites 
equilibrated with the bath occupation (fa). In Fig. [21 
this is manifested as the asymptotic normalized density 
(occupation probability) of 1. However, if pv = 0, tun- 
neling is forbidden and the system goes to an internal 
quasiequilibrium, i.e., the local excitation is spread in- 
side the system. In this case the asymptotic occupation 
probability of site 1 is 1/2. 



IV. STROBOSCOPIC REPRESENTATION OF 
THE INTERACTION PROCESSES 

Equation. (|37|) has two main difficulties for a numerical 
implementation: The first is the evaluation of the system 



8 



1.0 



0.8 



... 0.6 
9 0.4 
0.2 
0.0 




middle line 



bottom line p = 



v 



th 



8 



10 



si-: 



FIG. 2: (Color online) Occupation probability iGf ± (t) /h to 
find at site 1 a particle initially at site 2. Each line corre- 
sponds to different kinds, pv, of SE interactions. The plots 
correspond to x = VuTs^/h = 10 belonging to the swapping 
phase and f3 = 1. 



nonunitary propagators under inhomogeneous perturba- 
tions; the second is to keep track of all the previous states 
of the system to enable the integration over past times. 
We will show that the decay homogenization enables the 
implementation of an efficient numerical algorithm. First 
of all, we identify in expression (|44|) that e~ t ^ 2TSE ^ = s (t) 
is the system's survival probability to the environment in- 
terruption, i.e., the probability that the system remains 
coherent, and di;/(2rsE) = q (ti) dt; is the "interruption" 
probability in a differential time around ij. The inter- 
action of the environment is discretized in intervals r s t r . 
where it acts instantaneously. This stroboscopic interac- 
tion leads to 



where 



S (t) = (i- P r (t) , 

oo 

q(t) = 22 p 6 (t - mT stx .) , 
n (t) = int (t/T s tr.) ■ 



(67) 
(68) 



(69) 



Here, the stroboscopic interruptions may occur at the 
discrete times mr s t r . with a probability p. At time t there 
were n (t) possible interruptions. In the joint limit r str . — > 
and p — > such that 



p/r str . = 1/ (2r SE ) : 



(70) 



we recover the continuous expression (see the Appendix) . 
Introducing Eqs. (|6~7| and (f6"8"|) into the reduced den- 



sity expression ((44)) we obtain 



G< (t, t) = ft G (t, 0) G< (0, 0) G (0, t) (1 - p) 



,n(t) 



• / d£iT SE 5 (tj - t m ) 
J o m=i 



x G 0R (t, U) S< m . ( ti ) G 0A (t u t) p (1 - P ) n{t - U) 
and rewriting we have 



(71) 



G< (i, t) = h 2 G m (t, 0) G< (0, 0) G 0A (0, t) (1 - P ) n 

n 

+ h 2 Y / G 0R (t, t m ) (JG^. (t m , t m ) G 0A (t, t m ) 

xp(l- P r m , (72) 



m — 1 



where 



i(t), t m = mr str ., and 



t2 ^sym. (^) ' 



SGf n ,(t,t) 



(73) 



In this picture, the evolution between interruptions is 
governed by the system's propagators 



G OR (i,0) = --exp[-iH s t/n] (74) 



and 



G ,,A (0, t) = G m (t, 0) f 



(75) 



The spin bath stroboscopically interrupts the system evo- 
lution producing the decay of the coherent beam. This 
decay is compensated through the reinjection of proba- 
bility (or eventually of coherences) expressed in the in- 
stantaneous interruption function SG^ (t, i) which also 
contains actual injection (decay) from (to) the bath. 

The first term in the rhs of Eq. (|T2[) is the coher- 
ent system evolution weighted by its survival probability 
(1 — p) n . This is the upper branch in Fig. [3J The second 
term is the incoherent evolution involving all the deco- 
herent branches. The w-th term in the sum represents 
the evolution that had its last interruption at mr s t r . and 
since then survived coherently until iiT str . . Each of these 
terms is represented in Fig. [3] by all the branches with 
an interrupted state (gray dot, red online) at the hier- 
archy level m after which they survive without further 
interruptions until nT str . ■ This representation has an im- 
mediate resemblance to that introduced by Pascazio and 
Namiki to justify the dynamical Zeno effect [30| . 

As mentioned above, the solutions of Eqs. (|T2"|) and 
pi)) are both computationally demanding since they in- 
volve the storage of all the previous steps and reiter- 
ated summations. Thus, taking advantage of the self- 
similarity of the hierarchy levels in the interaction with 
the environment, we rearrange expression (|72[) into a 
form optimized for numerical computation, 



¥2'G < (t n +i,t n+ i) — 



G 0R (t n+1 ,t n ) G< (t n , t n ) G UA (t n , t n+1 ) (1 - p) 



OA 



G 0R (< n+1 , t n )6G< n] (t n ,t n )G UA (t n ,t n+1 )p. (76) 



OA 
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FIG. 3: (Color online) Quantum branching sequence for the 
stroboscopic evolution. Gray (red) dots represent states with 
interrupted (incoherent) evolution while the black dots are co- 
herent with their predecessor. The horizontal continuous ar- 
rows represent the isolated evolution and the vertical dashed 
lines are the instantaneous interruptions. Notice the self- 
similar structure. 



the system Hamiltonian Hs is 

h s = fi^ L (if + U) + \b 12 (/+f 2 - + 

= ^ (l{ + I*) + 6 12 (ff/f + PJ y 2 ) . (77) 

Here, the first term is the Zeeman energy and the second 
term gives a flip-flop or XY spin-spin interaction. The 
environment Hamiltonian is described by 



i>3 



i>3 



and for the SE interaction we have 



ftsE - a 23 I 2 I 3 + 5&23 ( 



(78) 



(79) 



where this spin-spin interaction is Ising if 623/023 = 0, 
and XY, isotropic (Heisenberg), or the truncated dipolar 
(secular) if 023/623 = 0, 1, —2, respectively. 

We map the spin system into a fcrmionic system using 
the Jordan- Wigner transformation (JWT) [33j |. 

f+ = c+exp|i 7 r^c+cA. (80) 

The previous Hamiltonians become 

= fo^L (ci"Ci + c+c 2 - 1) + |6 12 (cf c 2 + c^c^) , 



This equation provides a new computational procedure 
that only requires the storage of the density function at 
a single previous step. Besides, it avoids random aver- 
ages required in models that include decoherence through 
stochastic or kicked-like perturbations [U This 
strategy is being implemented in our group in various 
cases involving quantum dynamics of many spin systems 
in the presence of dissipation processes and decoherence. 
Equation (|7^)) manifests that the fictitious self-energy, 
proportional to (i n , t n ), compensates instantaneously 
the fictitious leaks and injection of Eq. (|40|) . The concep- 
tual consistency of the approach is illustrated by choos- 
ing 5G^- (t n ,t n ) = G < (i„,i„): one recovers a coherent 
isolated evolution. 



V. APPLICATION TO SPIN SYSTEMS 

We apply this procedure to the spin system of Ref. [TtJ 
providing a first principle derivation of the phenomeno- 
logical equations employed there. We consider a system 
with M — 2 spins 1/2 coupled to a spin environment with 
the following Hamiltonian H = Hs + He + Hse, where 



(81) 

He = (c+a, - |1) + E hhj + > 

i>3 ^ 3 

j>i 

(82) 

^se = a 23 (c+c 2 - |) (c^c 3 - §) + 5633 (c+c 3 + ct c 2 ) • 

(83) 

Here, the system interacts with the environment through 
site 3 (the surface site of the bath). In the last Hamil- 
tonians, the terms proportional to the identity do not 
contribute to the dynamics because they only change 
the total energy by a constant number. This Hamilto- 
nian describes a standard crosspolarization experiment 
(swapping gate) in NMR pji In this experiment, 

site 1 is a 13 C and site 2 is a 1 H while the environment 
is a X H spin bath. The typical experimental Hartmann- 
Hahn condition (l7l . [l9| equals the values of the effec- 
tive energies at the 13 C and the 1 H sites to optimize the 
polarization transfer. The SE interaction has terms lin- 
ear in the number operators c 2 c 2 and c 3 c 3 , that only 
change the energy of sites 2 and 3, respectively. Thus, 
the Hartmann-Hahn implementation, compensates the 
change of energy produced by the environment through 
these linear terms. Finally, we have Hamiltonians equiv- 
alent to those in Eqs. |T][3]) where the site energies are 
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equal, and V%2 



u. 



(dir.) 



2:i 



and 



U. 



(exch.) 
23 



0. 



The spin dynamics of the system is described by the 
spin correlation function [2(| H3] as follows: 



Pi2(t) 



(tt oq .|/f(f)J|(0)|*E oq .) 

(*eq.|/|(0)/|(0)|* eq .) 



(84) 



which gives the local polarization at time t on the i-th 
spin with an initial local excitation on the 2-nd spin at 
time t = 0. Here, l^cq.) is the thermodynamical many- 
body equilibrium state and 



If (i) = e ™/hjz e -mt/n 



(85) 



are the spin operators in the Heisenberg representation. 
After the JWT, the initial local excitation on site 2 is 
described by the nonequilibrium state 



(86) 



In the experimental high temperature regime, k^T much 
larger than any energy scale of the system, the spin cor- 
relation function becomes 



P i2 (t) = fG5(M)-l. 



(87) 



Notice that G^ (t, t) implicitly depends on the initial lo- 
cal excitation at site 2. Here, (t, t) is the reduced 
density function of sites 1 and 2 and can be split into 
the contributions G u (JV) (^2, *i) from each subspace with 
N particles (or equivalently N spins up) in the following 
way [27l ]: 



M 



N=l 



(M 

2 M-1 



G 



< (W) 



and analogous for the hole density function. The initial 
condition in this picture is described by 



G< (N >(0,0) 



N-l 



M- 



M-N x X 



(89) 



Finally, the solution of the polarization Pi2{t) is the same 
as that obtained in Eq. ([55| . 

By substituting in the present microscopic model 
Txy *-* Ty and Tzz IVi we obtain the same dynam- 



ics as that found in Ref . [17j for a phenomcnological spin 
model. There, we showed that such a solution presents 
a quantum dynamical phase transition in fair a gree ment 
with the phenomenon observed experimentally [18j . 



VI. CONCLUSIONS 

We have shown a method that involves the transforma- 
tion of the density function expressed in the Danielewicz 
integral form into a generalized Landauer-Buttiker equa- 
tion. This was possible by resorting to Wigner energy- 
time variables to perform the fast fluctuation approxi- 
mation for the environment which leads to interactions 
local in time. Further on, we effectively symmetrized 
the system-environment interactions transforming them 
into a spatially homogeneous process. This has a uni- 
form system-environment interaction rate leading to a 
simple non-Hermitian propagator. The original multiex- 
ponential decay processes are recovered by an injection 
density function. Moreover, through discretization of the 
GLBE, we built a stroboscopic process which is the basis 
for an optimal numerical algorithm where the quantum 
dynamics is calculated in discrete time steps. Finally, we 
applied these techniques to a spin system giving a mi- 
croscopic derivation that justifies the stroboscopic model 
used in Ref. [ItJ to explain the experimentally observed 
quantum dynamical phase transition. 



APPENDIX A: RECOVERING THE 
CONTINUOUS PROCESS 



In order to recover the continuous expression (|44p from 
the stroboscopic one ([72")) we notice that if n (t) = n, we 
can write Eq. (|6"T|) as 



where the first term is the equilibrium density (identical 
occupation for all the sites) and the second term is the 
nonequilibrium contribution where only site 2 is excited. 
Thus, we have an expression such as (fl"2)) for each iV-th 
subspace (see Ref. |27j ) . For this two-spin system, as we 
showed in [27j , the — 1 term of Eq. ([57)1 is canceled out by 
the background evolution, i.e., the evolution of the first 
term of Eq. ([89| plus the evolution of the second term 
of Eq. ()12[) for the N = 2 subspace. As a consequence, 
the observable dynamics only depends on the initial local 
excitation at site 2, 



G< (1, (0,0) 
and evolves in the first particle subspace, 



P i2 {t) 



h r < (i) 



(t,t). 



(90) 



(91) 



a(i) = (l-p) ™ 
If t = nr s tr. then 



\ (nr sta . )/T 8 tt. 



»(«)= ! 



r str. 



t/Tstr. 



(A2) 



By taking the limit r str . — * the variable t becomes con- 
tinuous yielding 



s(t)= lim I 1 - 

r atr .^0 V 2r SE 



t/T B tr. 



= exp [-*/ (2r SE )] , 
recovering the continuous expression for s (t) 



(A3) 
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By substituting p = T str ./(2rsE) in Eq. (|68|) we have The continuous expression of the GLBE (|44|) is then ob- 
tained. 

1 °° 

q(t) = V T stI .S(t - mr st r.). (A4) 

m— 1 

In the limit r str . — * 0, t m — TOT str . becomes a continuous 
variable and we can convert the sum into an integral, 
leading to 

1 f°° Hf 1 

q(t) = — / r stl 5(t ~ U— = (A5) 

<""SE 7o r str. ^SE 
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